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An Effect of the Atmospheric Electric Field upon 
the Observation of Ions in the Atmosphere 


By G. ISHIKAWA, M. KADENA and M. KITAMURA 


Meteorological Research Institute, Tokyo 


Abstract 


In the observation of ions in the atmosphere, we must consider 
that the observed values should differ from those in the free atmosphere 
owing to the deformation of the atmospheric electric field caused by the 
existence of the building, within which the observing instruments are Settied. 
Therefore, comparison among the observed values at various positions 
would be meaningless, without taking account of such an effect of the 

' electric field . 

In this paper, we tried to find the suitable position of the aperture 
for drawing the air into the building in a few simple cases, so that the 
consideration of such an effect may be dispensed with. — oi LAs 

1. The case of plane sien a iid ee.” : SN 
2 First we treat of the stationary state of i jonization on the plane ground, hich we 
consider as the standard case. Assuming that a uniform electric field is directed 
downward vertically, that only small ions exist, that the rate of ion production qis 
constant, and that the loss of ions is only due to the ae process, aes Te 


wing. 3 ee 


si mptions do not bsiiptantially sees our cde ane so ‘the ‘follow ing. expressio 
be. 


Ae ee epee te = 
; oe ilies So ae 2a 
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q—an*n- + hE SE =0 (1.4) 
q—on-n*— KES =0 Mess (1.5) 


from these equations, we easily obtain 
n*-+-n- =n=constant (independent of z and kE) (1.6) 


from (1.4) and (1.6), making EO 
; rh; Hig 
n*(n—n*) =~" and 2*=n" =) (1.7) 
SO, hae Ge =/q/a 


a te a into (1.4), we Lae 


(nt —/ aja qiap +E oh” <0 (8) 
from. this equation, - ] 
: nt= Yq/at+ oar (1.9) q 
Similarly, 

n- =V Wl aae0 (1.94 


where c is ‘the Beniteary integration con- 3 
_ Stant. F 
Taking account of tins. con- — 
ditions (#7),-0=0, (#*).-0=”, ¢€ is easily — 
determined, with which the ee 
rie ars (1.9) are written as follows, 


og) ee ip — : 
wn +3, gaa Slete) (1.10) 
= nate 
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E-yT=-1 (2.6) 
If; E.=E,=0, E,=—E at infinity, 
ne —E(1+ 4 sin 0 (2.7) 
= 
Ey=—E(1— =) cos 6 (2.8) 


where, F; and a are the 7 component of the electric field and the radius of the sphere 
respectively. 
From (2.6), (2.7) and (2.8), we have 
io Meee eater ss Pe 
(77+ a") sin@ ~ (#—a")cosé 


the solution for (2.9) is 


d0=EdT (2.9) 


ET(r,0)= (7c, C.) dr +ET(b,6:) (2.10) 
6 
where, 
CANE 4 be eee | 
fire.) =! (1+ =e oA (2.11) 
Lime | 
a= ee cos 0; (212) 


5 and 6, are the values of 7 and @ at the boundary, and the value of @ on the left 
side of the equation (2.10) is to be determined by the following equation after carrying 


out the integration, 


f (~P+a’) | 

6=cos JOT) cond Belo 

lbet+ar) n'y en) 
Since the elliptic integral flr, c,)dr is a little troublesome, we have integrated 

numerically. The result is shown in Fig. 2, 

where the following values are used, 


a=8nek (assume kt=k-) 
n=2X10*/cc 


cm/sec. 
volt/cm 


E=\1volt/cm 


= 


We shall next get the position Fig. 2 Equi-ion density line around a two 

at which the distribution of ions is dimensional semi-sphere. 

similar to those near the plane ground. There may be the same distribution of ions 
as that on the plane ground near a position on the boundary of the building, on which 
the same electric field as that on the plane ground exists. In order to get such a 
position, we have only to find the position on the boundary at which the relation of 
| E'| =E is satisfied, considering that the size of the building is large, and that the 
effect of the wind is now neglected. The substitution of this relation gives the. follow- 
ing result, aS 

2Esind=E or  4=30° 
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Therefore, in the case of a semi-circular cylinder-shaped building, the effect of the 
deformation of the electric field due to the building need not be taken into account, 
if the aperture is set at the position whose angle of elevation in 30 degrees. 


3. Two-dimensional elliptic boundary 

3. A Transversal case 

Since the three cases to be stated below are based on the same idea as the case 
of two-dimensional sphere, we shall omit detailed descriptions. We take the elliptic 
co-ordinate system in which the surfaces y=constant and =constant show the equi- 
elliptic surface and equi-hyperbolic surface respectively. Taking the surface, y=», as 
boundary and Ev, Ex as the v and » component of the electric field respectively, 
then 


1 adi eee j 
ye Se Sek 
“2 cosh*y—cos*p or es 
E,=— ee COS yt (B82) ; 


21 cosh*v—cos?y 


here, v and yw are connected with: the rectangular co-orninates as follows 
x=2acoshy cos p, “y=2asinhy sin pz (3.3) 
The position thus obtained on the boundary at which | | =E is satisfied is as follows 
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for example, taking Y/X=2/3 


0o==44°50/ 
4. The case of three-dimensional spherical boundary 


With the usual polar co-ordinate system, we write the electric field as follows ; 


3 
E,= Ecos0+2E, cos0 (4.1) 


4 24. So 
E,=—Esind +E, sin 0 (4.2) 
By the same process as in previous subsections, we easily find, 
Sasn= 1 ‘ OnO/ 
#,=sin“} 3 =19925 


5. Summary 


The suitable observing position, 0, dispensirig with the correction, as shown above, 
is aproximately determined by the geometric shape of the building and independently 
of k and E. In the case of a square shaped boundary, the calculation of E is difficult, 
but nevertheless, if we set the aperture for drawing the air on the position, where 
| F' | =E is satisfied, the observed value needs no correction, by the same consideration 
as above. 

Though we have discussed so far the case of no wind, we shall, finally, add the 
consideration of the effect of the wind. 

In this case, the expression of V* and V~ becomes as follows 

Vt=kE+V, 
V-=-kE+V, 
where VY, is the wind velocity. And assume that div V,=0 is Satisfied for simplicity. 
From (1.1), (1.2) and at stationary State, 
div (2*V+—n-V-)=0 
from this, kE-p (n*+n-)=Vo-7 (u*—n-) (5.2) 


(5.1) 


Therefore, the relation of m+-+~=constant is not always satisfied, and So it is clear 
that the distribution of ions in this case differs from that in the case of no wind, but 
in the windy case, to get an exact solution is not easy, therefore in place of doing 
so, we shall show the limit of the wind velocity within which the previous results 


may be used. Considering the effect of the wind, 
0 On- 
div n-V-=—kE-pn-+Vo-p-n- =—kE, oe + Voz 
y and tc show the direction of normal and parallel to the boundary respectively, and 
E,, Vo are the perpendicular and parallel components to the boundary of # and V oe 
From the previous results, only ions near the boundary contribute to the observed 
values. For example, we consider a two-dimensional sphere. At the boundary, the 


above equation is written as follows, 


ar V, 0 
div n-V-=—hE, 2 = Tae hy 63) 
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In this expression, if the first term is far larger than the second the effect of the wind 
may be neglected, 
namely, ay? Ri . F 
gob 65.4) 


On- 
| hE | aes 


re 
It must be noted that the value at the boundary must not be used in calculation, since 
the observed ions are not only those on the boundary but rather their space average 
value near the boundary. The distance from the boundary determining the region 
for the average depends on the velocity of the inhaled air and the mobility of.ions. But, 
here, in order to know the state roughly, we take the approximate value of ~ at the 
unit distance from the boundary for the above inequality. 


Then we have, 
| Vo | <0.4 acm/sec. 


where a is the value measured in meter. 

For example, if @=5 meters, the right side of the above is 2m/sec. So, at a 
wind velocity less than the order of 1m/sec., the results obtained in the previous 
sections may be approximately valid. Rigourous treatment in the case of strong wind 
will be performed in future. | 


__ . We should like to express our, hearty thanks to Dr. Hatakeyama, Director of 
Macmiancl Research Institute, for his kind suggestion in this work, 
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A Theory of Diurnal Magnetic Variations in Equatorial 
Regions and Conductivity of the Ionosphere E Region 


By Motokazu HIRONO 


Geophysical Institute, Kyoto University 


Abstract 


In two previous papers, I suggested that the electrical conductivi- 
ty in E-W direction of the £ region in a narrow belt along the magnetic 
equator is much greater than in the other latitudes, on account of a 
polarization field produced by the Hall Current, provided that 2 is 
much less than 10. And that this high conductivity will cause enhanced 
diurnal magnetic variations near the magnetic equator. 

In this paper, I discuss this theory further in detail. A dynamo 
theory for the F region with anisotropic conductivity, which is linked 
by highly conductive lines of magnetic force to the F region, is 
examined. A case that the E layer only makes a’ tidal oscillation, 
which produces observed lunar diurnal current system, is calculated. 
It is found that, in middle latitudes, the induced polarization field in 
E-W direction by the Hall Current is greater than that by the direct 
current, and the vertical drift of the F2 layer by this polarization field 
seems to be adequate for the observed F2 lunar variation. According 
to the present theory, conductivity of the EF region at the magnetic 
equator decreases with increase of the main magnetic field, and the 
distribution of range of magnetic variations along the magnetic equator 
reported by Egedal, seems to support this relation. 

1. Introduction . 
In the course of statistical yg  wancayo Losop 


researches of Sq field of terrestrial sgh 
magnetism in middle and low Es 
latitudes, it is noticed that the +!*° ae 
range of diurnal magnetic varia- ; x | 
tions at Huancayo near the mag- aes = 
netic equator (Dip 2.°3N) is much 100 


greater than at other stations. 7 
_ Whether there are other such 


13 


a ee 
150°E.M.T. 

Fig. 1 Diurnal Variation of Magnetic Horizon- 

places than Huancayo, is a matter : tal Intensity. Mean Values: Losop Feb. 

‘ 12-15, Huancayo Feb. 11-14 (1934). 

. of great interest, For the purpose , o)y © (After S. Hayami and. H) Higashinakay.9) 0.7 
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of studying solar eclipse of Feb. 14, 1934, Prof. S. Hayami and H. Higashinaka [1] 
made magnetic observations at Losop Island (6°53 N, 152°44/ E, Dip 1°45’ N) from Feb. 
12 to Feb. 15, which were reported by Prof. M. Hasegawa [3]. In Figure 1 mean H 
variations are reproduced for Losop (Feb. 12-15) and Huancayo (Feb. 11-14). 

Mean 4H of these are 1337 (instantaneous) for Losop and 1387 (hourly mean) 
for Huancayo. M. Hasegawa and M. Ota [2], made detailed analysis of the Sq field 
for the II polar year data and concluded that a line of maximum diurnal magnetic 
variation (in place of which we write MS line in the following) is supposed, and lies 
nearly along the magnetic equator (dip zero line). A.T. Price [4] suggests, using the 
above mentioned data, that this line has some deviation from the magnetic equator. 
But according to Hasegawa’s opinion, it is dangerous to decide this deviation from 
the above data, because observations near the magnetic equator have been very few, 
and direct observations in this region will be able to decide an exact deviation. 

J. Egedal [5] showed the ranges for 1948-50 near the equator, his results are 
reproduced in Fig. 2 in which ratios of the ranges to that of Huancayo are shown for 
South America and India. 


India 
(1950) 


South America 1.0 
(1949) 


Q8 08 
wl) 
3 2 
za < 
= 06 4 06 
AG 
04 04 
20°S [SSEMAG Oe 5° O° 5° MG ISN 
GEOGR. LATITUDE GEOGR.LATITUDE 


Fig. 2 Comparative Ranges of H refered to Huancayo (After A.A. 
Giesecke and J. Egedal) 

M: Magnetic Equator, G: Geomagnetic Equator (dipole equator) 
From these results it is supposed that MS line in the Far East and in the Pacific 
seems to lie very closely on the magnetic equator. Mean ranges of H at the magnetic 
equator shown by Egedal are as follows: South America 116;, India 105y, Togoland 
1247. These values are about as twice as that on the equator obtained by harmonic 
analysis. Magnitudes of main magnetic field are 0.37, 0.387, 0.3297 respectively for 
these places. South America and Togoland can hardly be compared because the 
latter’s only is of instantaneous value. The range at India seems to be smaller than | 
at South America and the main magnetic field of the former is greater than the 


latter’s. At present it seems promising that MS line agrees with the magnetic equator, 
therefore we examine this fact theoretically. 


2. Influence of Hall Current 

In two previous papers [6] ‘and [7], hereafter refered to as I and II respectively, E 
we made preliminary calculations on the influence of ‘the Hall Current (in place of ‘ 
which ‘we write H.C, in the following) to the electrical conductivity of the ionosphere, 
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The present paper describes a further examination of the problem. 

To discuss conductivity in the’ ionosphere, it is necessary to consider a mixture 
of four gases—neutral molecules, positive ions, negative ions, and electrons. These 
are distinguished respectively by suffixes a, +, — and e, and ions are denoted by 2. 
The mass, charge (in e.m.u.), velocity, number density of 7 th molecules are denoted 
by m,., @,, C,, 2,, and charge of an electron is specially denoted by —e. The number 


density, temperature of the gas as a whole are z and 7. The mass velocity ¢) is 


defined by Co = n,m, c,./ S\2,.M,, a 


and the peculiar velocity C.. is defined by C.=c,—¢e. The gas as a whole is supposed 
to be almost electrically neutral; 
thus >} 7,é,~0. (2) 


The electric current density J can be expressed in the form 


J= 3) n,e,€, = 330,€,(€o+C,) = y11,e,C, (3) 


D.R. Bates [8] suggests that 2_/nz,=A<1 in the EF and F region. But this value in 
the E region leaves some doubt, because reactions concerning molecular oxygen is 
not yet clear. T.G. Cowling [9], discussed the absorption of radio waves and suggested 
that 4<50 inthe F region. According to Cowling, when 4 is not much greater than 
unity or when 2,/2,>5.%10", the diffusion velocity C. can be expressed by the same 
form as for the free path method. In the ionosphere, the above conditions will be 
satisfied. When we discuss Sq current system, according to I § 1, we get the next 


relation as a good approximation, 


C,=(e,/m,){v,/(v,2+0,2) JE. — (1/402 (v,2+ o,”)}R x B+ (e,/m,v,)E 7 (v=-a) 

C,~0 (4) 
where E=E™ +c, x H (in e.m.u.), H“-=electrostatic field, H=the earth’s main mag- 
netic field, A=unit vector in the direction of H, 1/v,=m,D,/kT, D,=Diffusion coeffici- 
ent between 7 th molecule and neutral molecule. o,=e,4/m,, E., H,=Component 
vectors of E perpendicular and parallel to the magnetic field. using (2), (3) and (4) 
we get b : 
J=0B.+07 AX E+0,-Ej, . . (5) 
where aT = OVI, / my) Yn/ (Ye + wy) } ; 

a =(n,e/H){o."/(v2 +o") — 0.7/ (vi? +04") } 
Gn =€° S\Mr/Mr¥y) 
According to the dynamo theory, both horizontal tidal velocity ce, and the earth’s 
magnetic field H in the ionosphere, produce primarily. dynamo field c) x H, the vertical 
component of which is canceled by the resulting vertical electric field. The horizontal 


component of which produces a diurnal current system, together with a polarization 


field due to the resulting charge distribution. Near the E region, when 4<1 is satisfied, 
it is shown that o!~3x1072!2,, o%~5x10-%n,; therefore o'’/o'~10. Electric field 
produces not only a direct current, but also. considerable H.C, which is shown by the 
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second term of (5), therefore the dynamo theory may have to contain o” and differ from 
those of A. Schuster and S. Chapman. When 4 is not smaller than 1, ¢/o' becomes 
about 1/(24+1) times the above mentioned value, since a! becomes about (2A+1) times 
and o” is not increased. When 41 is satisfied, we can neglect H.C., and even when 
2 is small o/s’ rapidly decreases upward from the E region. At first we treat the 
case that the axes of rotation and magnetism of the earth coincide. 

Suppose that a uniform electric field is impressed in E-W direction near the equa- 
tor, then H.C. tends to flow vertically, but it cannot flow counteracted by a resulting 


polarization field; and as a resultant conductivity in E-W direction, we get 


aT = 9"{1+(a"/a")?}, (6) 
which is shown in II, (10). 

In this case positive and negative charges can move vertically by the same 
amount in the same direction. Let rectangular axis O-xyz be taken, so placed that 
ox, oy, and oz are directed to the south, east, and upward respectively and unit 
vectors in the direction of ox, oy, and oz be i, j and kK respectively. 

We treat first an ionosphere near the equator. Suppose that the layer is 
initially confined between z=a, and z=a, (a:><a,), electron density depends on z only 
and is mainly controled by electron production and recombination, and that A is 
uniform, gas density distribution is the same as the real EF region. This treatment 
will be adequate, since the thickness of the E region is several ten kms and the 
diurnal electric field is supposed to be approximately uniform over several thousand kms. 
The charge density is denoted by p, the electric field is supposed to be perpendicular 
to the magnetic field, then we have the next equation in this layer 


J=0'-E+oUhxE, (7) 
y-E=4ne'o, (8) 
0p/0t=—p-J (9) 


Suppose that constant electric field E) is impressed in oy direction from the time ¢=0, 
then H.C. produces polarization field E.k, and we get E=E.j+E-.k. Substituting this 
for E in (7) we have Jy=a'-E, to" -E,, J,=0'-E,—a" - Ep. (10) 

As electron production and recombination are so rapid that variations of o! and o”/ 
with current are negligible, the next equation is obtained 0J,/0t=0'-(0E./0t), on the 
other hand 0F,/0t=—4nc*J,, therefore 0]./0t=—4nc’o'J(a,<z<a) .°, Jz=Jne!*, a= 
1/47mc’s’. As in the E region we have o’~3x10-%!n, for gas density n=10" and A4<1, 
therefore a~3x10-*/n,. For so low an electron density 2,=1, follows a~3x10-?sec. 
For n,=10° we have a=3x10-'sec, but it is doubtful that we can treat the processes 
in such a short time by our present theory, since 1/y;~2x10-‘sec. at m=10'. Never- 
theless we can deduce that J, rapidly tends to zero after a time, say 2x 107*sec. When 
A is not smaller than 1, a becomes about 1/(24+1) times the above value. Thus almost 
instantaneously, a steady state is attained and the relation J,=o'E,—o"E)=0 (1) 

holds at any height of the layer. At the same time from (10) J,=0" {1+(0""/a")"}-E,, 
therefore the conductivity of equation (6) is obtained at any height. Using equation 


’ 
< 
- 
cS 
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(7) we get 7-J=o'7-E+o"-p(hxE)+po'-E+po"-hx E. Writing o'=n,-0,', o@= 
n,-a,", and considering equ. (8) and V-(hx E)=—0E,/0z2=0 we get the following 
equation 

—00/0t=4rc?-a'- 9+ n4E(00,' /0z)—E,(00,1 /0z)}+(o,'-E,—o0,!' «Ey)-0n,/02z (12) 
At the steady state, it follows from this equation that, 


0 
Oz 


o=(E . a oy —E, ay!) Arco! (13) 


This quantity depends on gas density and not on the electron density. In the F region 
we infer that oy’ ~(24+1)-1.76 x 10-8/n, 

0,74 ~(1/H)1.59 x 10-7°(1—1.12 x 10° /n°). 
We take local scale height as 10km, then 0n/0z=10-°-n. We have from (13), using 
these relation, p= —E,-9.0310-*8 for 2<1, and n=10'%. When solar tidal oscillations 
in this layer are /, times as great as those at the ground, we get E,)~3.7/, (e.m.u.) on 
the equator according to a later section, then we get 

An, =/(—e)=2.17 x 10-7-1,. 
As 7, is estimated to be the order of 10°, 42, is very small. Thus we can see how 
well the condition of electrical neutrality (2) is satisfied. When A<{1 holds, we can 
also see the the same result. Next denote the total charge in a vertical. column, 
which has a unit cross section on a plane z=z, and is confined in the range z<2:, by 
q, then we have at z=z, E.=4zc*g. At the steady state, using (11) we obtain q= 
(1/42c*)(a" /o')E,=5.8x 10-7! 2, for n=10'3 at z=z,. Let the length of the column be d, 
then from the above relation it follows that d-4n,=q/e=3.65 x 10 for 1,=10°. For the 
value d~10km, we get 4dn,~3.65 x 1075. ; : 

If gas density in this layer is uniform, space charge distribution p in the layer 
vanishes almost instantaneously and EF, arises by a surface charge distribution. q. 
From (11), as the mean value in the E region, we get E,=6-E) for 4<1. 

3. Conductivity at any Latitudes 

We treat the same the distribution of ionosphere as in the previous section. 
Here the magnetic field is supposed to be uniform with dip ¢ (constant). When a 
uniform electric field is impressed in either E-W direction or in N-S direction, we can 
show that the vertical current vanishes almost instantaneously on account of a vertical 
polarization field, but in this steady state, a vertical mass movement of the layer 
remains. When a uniform electric field (E,%+E,j) is impressed, in the steady state 


we have 
bs =6,F +62 Ey, Jy =Cye kt oyky ; (14) 
E,={o" cos'¢-E,—(¢,—a')sin ¢cos ¢-E,}/(o,-sin’¢ + 07+cos’¢) (15) 
where g,=0'/{sin?¢+(a7/o ,)cos*} . 
6,y=0''sin ¢/{sin?¢+ (a"/o ,) cos’¢} (16) 
Oye = — Oxy — 


Byj= ot of (a cos ¢)?/(o7cos*$+<a ,Sin*¢) 


The values of «, and , were already discussed in I, § 3 but o,,y in I, (17) is not valid . 
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Fig. 3. The calculated values of 
(a) o,/o7, (b) o,./o with latitude, 
and (c) o//m,, otI/n,. with gas- 
density. 
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since it does not take the vertical polarization field into consideration. The values e 
of a,/s', ¢,/o' are illustrated with latitude in Figs 3 (a) and (b). o’/n, and o'’/n, with — 
gas density are shown in Fig. 3 (c). Here we use the value by Ferraro [10] for D: ‘ 
with molecular weight 20 for ion and 25 for neutral molecule. D, is refered to ga ; 
kinetic value, and AH is taken to be 0. 3 (em.u.) Therefore, we take the values goa) 
131x10°-n, v,=4.36 x 10-"-n. At middle latitudes o, a. and Oy~o". The a 
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An ideal case that these two layers are electromagnetically independent, was 
first raised by Dr. K. Maeda Ldide}s | 

But we think that interactions of these layers are important. At present it is 
most probable that the Sq current flows mainly in the F region. We imagine the 
following model for Sq field in the daytime. In the northern hemisphere, anticlock- 
wise electric field (dynamo field plus electrostatic field) of a whirl shape will be 
present. Direct curect will flow along this field as shown by I. Fig. 4. Horizontal 
component of H.C. will flow in such a way as to concentrate negative charges to the 
focus of the whirl. This polarization will, however, be considerably dissipated by 
direct current in the F region in the opposite direction to the lower horizontal H.C.: 
since in higher latitudes than about 7.°5, as shown in II, the E region is connected 
by highly conductive lines of magnetic force to the F region. Therefore the effective 
conductivity of the E region in these latitudes will be much less than that at the 
magnetic equator. 

In the following we suppose that 4<1, EF and Fi layers. are thin spherical shells 
with the above mentioned conductivities, and their distances are so small that they 
can be represented by a common radius a in mathematical treatment. 

Let & and Fi layers be Chapman layers with local scale height 10km and 30km, 
gas densities at these centers be 6x10" and 10", and maximum electron densities 
10°/ec and 2.5x10°/cc. For this model, we obtain a distribution of o, with height. 
Between FE and Fi layers we obtain by a numerical calculation, 1/(1/o,)=7.95 x 10-8 
(e.m.u.) which is very great compared with the values at respective centers, o(E)= 
4.4x10- and o(Fi)=3.58x10-". Therefore we can show that E and Fl regions are 
connected by highly conductive lines of force in middle latitudes. But vertical current 
will be very small. Here we assume that Sq current system is nearly horizontal, 
ie. J.<VJ2+J,. Considering the relation ¢,>0', 6%, equation (15) can be reduced to 

E,=—cot¢-£,. (19) 
On the other hand, as the third term of equation (5) can be written as 
o,H,=0—E,cos¢—E,sing)h 

equation (19) implies that #, is very small. Therefore, if B, deviates a little from 
the value of (15), current flows almost in the direction of lines of magnetic field, and 
the above assumption is violated. Thus equation (15) and consequently (14),-(16), (17) 
and (18) hold with good approximation, in this case. As the relation between F,, £,, and 
E, in the equation (15) is linear, we can divide £, into two parts; the one produced by 
the horizontal dynamo field, and the other produced by the horizontal electrostatic 
field. These sets satisfy the equation (15) respectively. In the following we treat the 
latter set. Leta vertical square abcd be taken, as shown in Fig. 4, so placed that ab 
and cd are horizontal and situated at the center of FE and Fl region, and in the 
direction of electric field in a practical sense, and that this square is fixed to the 
earth. Horizontal electrostatic fields on ab and cd be E, and £’,, which are measured 
positively in the direction of ab, upward vertical component of electrostatic fields on 


be and ad be E, and E’,. 
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Eds=0 around abcd, we obtain 


the following equation 
B,-E', = BY,—E,, (20) 


where bar denotes mean value on each side 
of the square. 

Let a length of one side be Z, which 
is of the order of 10’cm. Direction cosine 
Fig. 4 of ab be (J,m,0), then differential coefficient 


in the direction of ab can be written, using equation (19) as 

dE,/ds = —\(E.,-0 cot ¢/a00+ cot ¢-0E,,/a00)—m (cot 6-0E,,/asin 00) 
We deduce, from the observed results, that 

0E,,/a00~E,/a, OF ,/asin 00A~E,/a 
Therefore 0E./0s = —(E,/a)[1{2/sin’¢ (1+3 cos?0)+ cot ¢}-+m cot ¢] (21) 
Under the assumption that horizontal electric field in every region is of the same 
order of magnitude as that at the center of the region, and considering that E’,—E, 
~(dE./ds)L, we can make the following estimation. In higher latitudes than 20°, 


Se Pag i e 1 1 
[BB ISIE +1 Bil h 0 (t st ag) | 


|E-2i1<|B1+1/2n1- 


This relation holds approximately near 7.°5 but does not hold in lower latitudes if ab 
is in the N-S direction. Therefore in higher latitudes than 7.°5 we can take E,=E, 
as a first approximation. Thus we write Xi:=X:,=X, Yi=Y2=Y in (17) and (18). 
The above result will be approximately valid for E and F region. For such a 
uniform current system as Sq, difference of the height of current sheet, of 100km, 
gives negligible difference to the variation of terrestrial magnetism at the ground. 
Therefore, Sq is produced by J=4,+ 4 in a good approximation. We can, further, 
treat in such a way that 01,/a00+-0I,/asin004 =0, to be derived from 7 - J=0, considering 
J.=0 on the outer surfaces of the combined layer. 
Then the following two cases are examined. 
(I) F region is at rest ie. w,=v,=0 and E region makes tidal oscillations. 


We write J=/,i+1,j, .<I,, [,), and from (17) and (18), get 
a; i a,(uiH, +X)+4,,(—mH, =F Y)+8,X 

Ue 
[,, = @y2(1H, + X)+a,(—mH,+ Y)+8,Y ( : ) 


Xi+ Yj can be derived from a potential S, ie. X=—0S/a00, Y=—0OS/asin 002. 
Substituting these values in (22) and by the eiimination of S, we have 


yr 0/1 d /siné d /siné,\ 
ile ook) NA pa. | ae BALEY 
oe ) ale, “Iy)— o's, *)* aa i 1.) 


= BHM + M.v,)} spare og 1H: sin O(M3u;- aia Myv;)} f (23) 


et Sh 
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where 
K.= (a, +B, J. Ky = (y+ By) +9 
; Kay= @zy+(a,+8,)Gy+By)/aey,  Kye=—Koys 
Mi = —(By*G,y)/4, M, = {@,(ay+By)+A,y7}/d, 
Ms ={(4.+Bx)ty+Ozy7}/d, M;=(Bx*Gzy)/A, 
with A=(4,+8,)(Qy+By)+Ayy? 


Here we assume that electron density is constant at the same altitude. At the lati- 


tudes, here concerned, we get the following values approximately 
( K.=Ki/sin'g,  Ky=Ki=(a'+89+(a"8/(al-+B%) 
 Ky=Ki/sing,  Ky=a"?+(0! 4802/0! 


where a'=§,0'dh etc. (¢=Dip). 

M,=—7:/sing, M.=r2, Ms=r2,. Ms=7;/sin 4, 
where 11=B" a" /G, r2={a"(a'+B)+(a™)*}/G, 
with G=(a' +"? +(e") 
Ki, Ko, v1, 2 are independent of 4 and 0, if we neglect very slow variations depending 
on the magnitude of the magnetic field. Let velocity potential of the air be 
Vo=ks Pf sin (at+o%), (24) 
where kf is taken to be positive, ¢=i/+4 denotes local time, P” denotes associated 
Legendre function of Neumann. 
Then 1 =00%/a00, v,=0¢%/asin 002. 
We write H.=Ccos@ (C<0), and the right hand side of equation (23) ‘is reduced to 
the form, 


7 sin 7 ee ¢.1S8in (of +0%)+71-00(0)P% cos (oto)! (25) 
where 0 (0)=3 cos*0/(1+3 cos?0)-sin ¢ 


We shall treat, mainly the semi-diurnal air velocity (s=2). For numerical example, 
two cases (a), and (b) are considered. 
Case (a) oa =5K10~, a“ =3x10-*; pi=4K10- (ean) 
using these values we get 
| K,=13a!, Ky=19.50'; 71=0.41, 7.=0.385 
Case (b) ai =5x 107 ae 2107, Balox l0-? 
using these values we get 
Ki =13.60', Ko=8.16a"; 71=0.346, 72=0.77_ 
These values of case (a) are obtained in II. If Bates and Massey’s [17] estima- 
tion of gas density depending on rocket measurements, are correct, then case (b) will 


be better. At present the effective conductivity of the F layer to the electrostatic 


field is not yet clear and these examples will give rough order of magnitude to 
the polling. calculation. 


The values of tales and | 71 are 


16 M. HIRONO 


shown with colatitude 0, in Fig. 5. As 
the latter’s is much less than the former's, 
we can neglect the second term of (25), 
to discuss the order of magnitude. For 
case (b) we obtain the same result. 
Since p-L=0, I can be derived 
from a current function R, i.e. 
T,=0R/asin004,  I,=—OR/a00. 
After substitution, of these values, the 
left hand side of equ. (23) is reduced to 


the next form, 


1 aR @ /.» OR\| 
op sin 002° + 5g ( sing 00 ) f 


0 
OR |+ 1 2cosf OR. 


Hig *sind02 ak, "1430820 OA 
7 (26) 
Colatitude Suppose that the current function R, to 


Fig. 5 The Values of r2+(8/5)Ps°, 11-22(6)P2?;_ be determined from (73), is expressed by 
sin 0, cos*¢/sin 9, cos¢/(1+3cos"0); and t Lh . ical se = ; 

A) Beaks, coustion da pashvadate as {ek ge ee 

tude. (27) 

By substitution of this value, (26) is reduced to the form, 


RE sin 0-n(z+1)r% mP™ sin (mt-+o =) eee m? «7's Psin (mt + of} 


1 
aky 
cos ¢ 
— te *7aP% cos (mi+ ax) fin (28) 
As we are treating semi-diurnal variations, m=2, sa moreover ”=3, since the current. 
must be symmetrical against the equator. The values of sind, cos? *¢/sin 0, cos b/(L+ 
3.cos?@) are shown in Fig. 5 with colatiude. From these results we can eliminate 
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a OSA OR 
Mite “aod RK. asindod 
Pee oe bis alae (31) 
Wiis in OA OK “ao0” 


From these two equations, by eliminating S, and substituting (27) for R, and K, for 
K, the same equation as (29) is obtained. Hence K, is interpreted as an effective 
conductivity to the dynamo current in these latitudes. K, is represented as 
K,=a' +(a")*/(a' +B, (32) 
and K,=5a' for case (a), K,=10a! for case (b). 
(II) # and Fl regions make tidal oscillations with the same velocity i.e. 
ets, U2 = 01. 


In place of equation (22) we obtain 


'# = (a,+8,)(u1H-. +X)+ae,(—mH,+ Y) 
(33) 
L, = Gyx° (UH, + X)+(ay+By)(—aH,+ Y) 
From these equations we have, in place of (23), 
0 oa 0 1 0 (sind 0 /siné 
= ( ane ( : 
an xs") ~ a (ey) a0 (MK, *) + G0 SE") 
0 Chae 
ye (v,H.) +, (sin 0-1,H.) (34) 


The left hand side of this equation can be reduced in the same way as before, and 
the right hand side is reduced without any neglection, to the first term of (25), in 
which 1 is substituted for 7; By the same deduction as before, effective conductivity 
K, is obtained. The value of A; is Ki=K,;+£'. Therefore, in this case also, effective 
conductivity of the E layer is considered to be K,. 
4. Sq current near the Magnetic Equator 
In a narrow belt along line of zero dip of several degrees wide, the equation 

(19) may not hold, since both current and magnetic fields are nearly horizontal. The 
dy, in the equation (14) will not be valid and «a,, in I, (17) might be better and 

the difference of these values is very great. At present we have not been able to 

find the best value of o,,. Thus, we treat this region rather in the same way as in 

section 2. 

In this belt in the E region, we have 
TiewarXs, — yal’! -Y; 

where: a2 =(g0""dh, . as=Seo,dh. 

In this belt, in the Fl region, h,~=B;X., ty=B'Y2. 

By the same way as in middle latitudes, we can show that Yi~Y2, but we have 
Xi+-X,. Hence effective conductivity in E-W direction will be a’’7 + ' and that in 
N:S direction will be a certain average of a, and 8,. The latter will be larger than 

' the former, since a, is much greater than a’, o% ando”. As the most resistant part 
is effective to the dynamo current, a”“+ 7 will be the effective conductivity in this 


belt. 


“a — a. . 
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For case (a) all 4. pl =45a", (a *+89/Ki=3.5, a4 /K,=75 
and for case (b) alll 4 Bl =39.6a", (a !+89/Ki=3, a4 /K,=3.7. 
For any case, effective conductivity in this belt is greater than 3 times that in middle 
latitudes. This conductivity will slowly decrease from the equator towards 7.°S of 
latitudes. Therefore, to see the outline of the Sq current distribution, we may use 
the dynamo equation (31) for the following distribution of K, (A) A=; (const.) for 
0<0<0, (B) K=vKz for 0:<0<x—0, (C) K=K, for z—0,<0<2, 0,=85° is used in 
this section. From equation (31), by eliminating S we obtain 

0 (wH_sin @) } 


PR. 0 (5.9. 9R)_ axl WOH.) A 
~ sin 002s 08 (sing. a) aK 61 ohten sabes: CH (35) 
First, we consider the solution of the equation 
OR 0 OR\_ 
sin 002° +59 (sin a 00 )=0 (38) 


Substituting R=/(0)-F(4) in the above, we get equations (20) and (21) in I. The 
solution of (21) in II is-denoted by Qo” by ee [13]. Using his formula, an integral 
of (36) is expressed as 

SHAML—2) ee ee ee (37) 
where M is an integer, “.=cos 0, and A,,, By, Dy, Ey are integrating constants. A 
general value of the integral of (35) is obtained as the sum of a particular integral 
of (35) and the above integral. Thus, for semi-diurnal variation we obtain the follow- 
ing solution 
(A) R=B-KiAP3°(8)+Cigi(9)} sin 2é+-a) 
(B) R= BvKi{ P3*(P)+Crqi(9) + Csqe(@)} sin (2+ a) q 
(C) R= B-KiAP3*(0)+-Cig2(9)} sin 2t+a) reve.” we 
where B=(2/15)Ck,*, gi(@)=(1—#)/(1+2); q(O)= (1+ e)/—p) on Yerexs pel 
Considering the continuity of the normal component of, shietacrent sociale 
Se ieee Se ne Penh opoes erent hereon (AL ae 
wer igs sib.ctoie .\!) Sree bes bile ee stward 


AAIMS 


te ete 
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5. Vertical Drift Velocity of F2 Layer produced by the Hall Current 

In this section we discuss the case that only E region makes tidal oscillations. 
In section 3, we showed that E;~E,, and this result will be approximately valid for 
E and F region. Since this conclusion depends upon some assumption and deduction, 
it is clearly desirable to submit it to a direct test. It has been shown above that 
when we discuss semi-diurnal variation, in higher latitudes than 7.°5, the current func- 
tion R=rs;°Ps° sin (2é+a), a=a2" is a fairly good approximation for the current. 


From equations (22) and the above value we obtain 


Y=(1/a)-C-k,2{ A sin (2¢-+a)+B sin (2t+44-90°)} (38) 
=(1/a)-C-k,’- D sin (2t+a+¢), (39) 
where A=(K,/Ki){(—2/15)0P32/00+cos0-0P,°/00} 


B=(4/3)(K,/Ko) sin 6-P2'+(2/sin $)- P2'+(a7/Ki—a'/Ko) 


In F2 region the second term of equation (4) is the same for ions and electrons, and 
a vertical component W of this term represents the vertical drift of this region. Thus, 


we obtain 


Pix OS fay 
H.W 1+3 cos? 


(40) 


D.F. Martyn [15] suggested that main seat of the lunar current is the D region. 
Later, it was pointed out by Cowling and Borger [16], that conductivity of the D 
region is too low to be a main seat. The same objection was stressed by Bates and 
Massey [17]. Though the F region has considerable conductivity, induction drag [9] 
and kinematic viscosity are very effective in this region. Tidal motion in this region 
may be very small, from the theoretical point of view. Therefore, it will be possible 
that lunar main current flows in the EF region. Suppose that lunar main current flows 
in the E region, then a=—17° should be satisfied to agree with Chapman’s results for 
1905 equinox (14). 
If Hall Current is not taken into account, Y is of the same phase as the first - 
term of equation (38), then W is almost downward during 6 hours after lunar transit. 
While as regards world wide lunar F2 variation, except very near the magnetic 
equator, amplitudes of oscillation of h’F2 and #73* are about 2~4km, and maximum 
height is attained about 6 hours after lunar transit [15]. Therefore alike Martyn’s, 
a conclusion might arise that a reverse lunar currrent flows both in £ and F region. 
But for case (a), since | B| is greater than | A |, as shown in Fig. 5, and ¢ becomes 
about 70°, therefore upward maximum of W is attained about 7 hours after lunar 
transit. Hence no contradiction will take place with F2 variation, when the £ region 
is the main seat of the current. Pant 
When tidal oscillation in the E region is 1; times as great as that at the ground, 


amplitude of W at 0=60° is 0.527 Z,cm/sec. For case (b) a similar result is obtained. 


As the observations of lunar E variations are very few and detailed character-. 
istics are not yet clear, we have hesitated to discuss them. 
We tentatively suggest that the main seatof the lunar current is the £ region. 
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6. Variation of Conductivity depending on the Main Magnetic Field and 
We shall consider the variation of conductivity depending on main magnetic 
field and 4. Here we take the value of D; which was suggested by Tyndall (1938) [18] 


using observed mobility of N.* and that of D, which was shown by Huxley (1949) 


[19]. Molecular weights of ion and neutral molecule are taken to be 26, temperature 
of the gas to be 300°K. Therefore we take the values v,=8.4x 1079-2, v,;=5.09 x 10-M-n. 
We represent the E layer at the equator at noon by a Chapman distribution with 
center at the gas density 2=10" and for local scale height of 10%cm. 

According to examination in section 2, o’/‘ will vary to some extent depending 


7T varies considerably with 


on the magnitude of H, since o’ varies little with H but g 
H. We calculated conductivity of this layer using equ. (6), by numerical integration 
for various values of 4 and H. S; denotes the case H=0.3 and S; the case H=0.4 
(e.m.u.) 

On the other hand we calculate conductivity A of this layer using o7 and For- 
mula by Appleton [20] for the case A=0, and get the following values A=A,+Ay 
where A, =7,0-1.55x 10-"H * A,=N.o°3.1X10-"H-? with x. the maximum electron 
density of this layer. 

When A=+=0, we assumed that A= AORN As The ratio S/A with 4 is 


10 } 


Table : 
> cs | "10 | 20 
S/S) | 1.4 | 1.27 | “119 


illustrated in Fig. 7 for H=0.3 and 0.4. 
£2 The ratio S;/S; for various —— of « 

M=0.4 — ryt pike is shown in Table 1. oT 2 ae 
ica: The ratio S/A may be a little 

wy eth he. th Yet osmallerithan. curves in Fig: 7; pam 
: ee fields are net ompletely 


oa i Pt —_—_— 


"ne 
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conductivity in E-W direction due to the polarization field near the EK layer produced 
by Hall Current. The above calculation is made for the case that the axes of rotation 
and magnetism coincide. In the actual case that the two axes do not coincide and 
local anomalies are present, MS line will be along the magnetic equator (line of dip 
zero). If the concentration of the negative ions is effectively zero, the effective con- 
ductivity of the £ layer in middle latitudes is shown to be some 5 times that obtained 
by Appleton’s formula. It’ is tentatively suggested that the main seat of the lunar 
current is the & region, considering the vertical drift of the F2 layer, produced by 
the polarization field due to the Hall Current. According to the present theory con- 
ductivity of the E layer decreases with increase of the main magnetic field at the 
magnetic equator and observed distribution of ranges of diurnal variations reported 
by Egedal, seems to support this relation. 

Moreover this theory is effective only when 4 is much smaller than 10. Recent 
observed results seems to be favourable to the present theory. Further observations 
near the equators are very desirable. 


Acknowledgements 


The author wishes to express his thanks to Prof. M. Hasegawa for his continued 
interest and advice, and to Prof. T. Nagata for his valuable discussions during the 
progress of this work. 

Many thanks are due to Dr. K. Maeda for his criticisms and helpful advice in 
this work. 


References 


[1] S. Hayami and H. Higashinaka, J~ Shanghai Institute (in Japanese), Sec. I, 3, 
29 (1934). 
[2] M. Hasegawa and M. Ota, ‘Transactions of Oslo meeting, 1948,’ 426 (1950). 
[3] M. Hasegawa, ‘Transactions of Brussel meeting, 1951’ (in the Press). 
[4] A.T. Price, J. Geophys. Res., 56, No. 2, 259 (1951). 
[5] J. Egedal, ‘Transactions of Brussel meeting, 1951, (in the Press). 
[6] M. Hirono, J. Geomag. Geoel. Kyoto, 2, 1 (1950). 
[7] M. Hirono, J. Geomag. Geoel. Kyoto, 2, 113 (1950). 
[8] D.R. Bates and H.S.W. Massey, Proc. Roy. Soc. Lon. A, 187, 261 (1946). 
[9] T.G. Cowling, Proc. Roy. Soc. Lon. A, 183, 453 (1945). 
[10] V.C.A. Ferraro, Terr. Mag. Washington D.C., 50, 215 (1945). 
[11] K. Maeda, Read at the meeting of Ionosphere Research Committee, Feb. 29. 1952. 
[12] S. Chapman and J. Bartels, ‘Geomagnetism II,’ 768 (1940). 
[13] E.W. Hobson ‘The Theory of Spherical and Ellipsoidal Harmonics’ 89, Camb- 
ridge (1931). 
[14] S. Chapman, Phil. Tran. Roy. Soc. Lon., 218, 1 (1919). ~~ 
[15] D.F. Martyn, Proc. Roy. Soc. Lon., A, 194, 445 (1948). 
[16] T.G. Cowling and R. Borger, Nature, Lon. 162, 143 (1948). - 
[17] D.R. Bates and H.S.W. Massey, J. Atmos. and Terr. Phys.,, 2,1 (195). 
[18] A.M. Tyndall, Positive Ions, Camb. (1938). 
[19] L.G.H. Huxley, A.A. Zaazou, Proc. Roy. Soc. Lon., 199, 402 (1949). 
[20] E.V. Appleton, Proc, Roy. Soc. Lon. 162, 451 (1937). 
(Read Apr. 28, 1952) 


Self-Reversal of Thermo-Remanent Magnetism of Igneous Rocks 


By Takesi NAGATA, Seiya UYEDA and Syun-iti AKIMOTO 


(Geophysical Institute, Tokyo University) 
(Read Oct. 22, 1951 and Apr. 29, 1952) 


Abstract 


An experimental proof of self-reversal of thermo-remanent 
magnetization of igneous rocks is obtained. By cooling in a weak 
magnetic field, some igneous rocks and the ferromagnetic minerals 
separated from the mother rocks can have the remanent magnetization 
whose direction is opposite to that of the magnetic field during cooling. 
The processes of development and decay of the reverse thermo-rema- 
nent magnetization by cooling and heating respectively were directly 
observed (Figs. 11~19) and the characteristics of its appearance were 
examined. A few possible mechanisms of appearance of this peculiar 
phenomenon are discussed with reference to Néel’s theory. | : 

1. Introduction 

It has frequently been reported (1), (2), (3) that some igneous rock masses have 
the remanent magnetism, the direction of which is almost opposite to that of the 
present geomagnetic force. This phenomenon of reverse Senahent magnetism has 
been found at various localities in the world. 

On the other hand, that the remanent magnetism of igneous ee is copay 4 
due to the so-called thermo-remanent_ magnetism has been experimentally. proved by 
a number of investigators (4), (5), (6) and theoretically treated by L. Néel Ae i 
is, a mass of igneous rock after-cooling | ina weak magnetic field from a ter 


above the Curie point has generally a relatively strong remanent reciaieat Saar ah the 
direction of which is just eae with stat of f the ‘agnetiedeld erelicinag aia 
a; cooling. no accetglin se : grein te Tadtidies ao taptipics “OT: . wuts caeenrte a Se 
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However, here is a question: Js there any other process of producing the 
remanent magnetism of igneous rocks? Or more directly speaking, is there any 
process of producing the remanent magnetism whose direction is opposite to the 
applied magnetic field? 

L. Néel (10) has recently proposed a few theoretical possibilities of the self- 
reversal of thermo-remanent magnetism. The present writers have also been attacking 
the above-mentioned question, experimentally and theoretically. In this preliminary 
note, the results of experiment of demonstrating the reversal of thermo-remanent 
magnetism will be reported. 


2. General Description of the Test Sample 


The rock specimen concerned here is a hypersthene hornblende dacite of Mt. 
Haruna, an extinct volcano in Kwantéd District. These dacitic rocks are distributed 
widely on the north-east slope of Mt. Haruna, and they are of the form of pumice of 
fairly large size, 40cm in their mean diameter. 

These dacitic pumices have the natural remanent magnetism of 4.3 x 10-°¢.272.1. 
in average intensity, and its direction is almost opposite to that of the present geo- 
magnetic force, ie. D=215° [= —14°. 

For the purpose of examining in detail the magnetic characteristics refered to 
its chemical and mineralogical compositions, the ferromagnetic minerals contained in 
the .rock were separated by two different methods. i.e. by means of a magnetic 
separator and by a panning method. Those mineral grains are the phenocrysts in the 
mother rocks, being 0.4mm in mean diameter. The result of chemical analysis of the 


{x10 Table I ferromagnetic mineral by I. Iwa- 
Chemical Composition Saki and K. Katsura (11) is given 
; (AK 51070801 b) in Table I where the numerals 
FeO; | 53.15%] 57.994 in the righthand side mean the 
10) Reo | ates 347 relative amountes of Fe,O3;, FeO 
and TiO,, excluding the other 

8 TiO, 6.82 7.4 | nl asuanhee 
MnO | 0.69 sae The change of magnetic 
$ VzOs | 0.49 = susceptibility of the ferromag- 
, netic mineral in a weak magnetic 
field was measured by a ballistic method apparatus 
2 (12), the result being illustrated in Fig. 1. As will be 
seen in this figure, the change in susceptibility of the 
0 200400 coor. oT sample with temperature is stepwise, having appar- 
Fig. 1. Change of magntic — ently two Curie points, ie. 540°C and about 280°C. 
susceptibility with tempera- This result may indicate that the ferromagnetic 
4 i a ay pest minerals concerned here are composed mainly of two 


kinds of ferromagnetic substances, one being, of course, 


t magnetite, while the other might be a combined state of Fe:0;-FeQ-TiO,, | 
é 


i : 
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.& On the other hand, Fig. 2 shows the 
ou intensity of saturation magnetization of this 
mineral in the temperature range from 10°G 

to 650°C. The magnetization changes with 


tg temperature almost in accordance with Curie’s 

Law, where Curie point is about 570°C. It can 

40} be said, therefore, that the largest part of the 
ferromagnetic mineral is magnetite. 

3. Reversal of Thermo-Remanent Magnetism 

) 200 400 600% T ; 

Rigs 2. (Change,Ot .sattatonemaenets of Rock Sample and Separated Ferromagnetic 

zation with temperature. Minerals. 
(= 2000 Oe.) The procedure of magnetic field cooling 


of test sample were carried out under these different conditions: i.e. (a) in a non- 
inductive and non-magnetic vacuum furnace, (b) in a non-inductive and non-magnetic 
N,-furnace, where the pressure of N; gas was kept around 780mmHg, and (c) in a non- 
magnetic furnace in the free atmospheric air. The magnetic field in the furnace can 
be varied with the aid of a Helmholtz’s coil. 

The test samples of the rock were of the form of a circular cylinder, 10cm long and 
2cm in diameter, while the ferromagnetic grains were put in a cylindrical tube of 
fuzed silica, 10cm long and 6em in its yu 
inner diameter. 

The thermo-remanent magnetiza- 
tions (T.R.M.s) produced by cooling 
down the test sample to the room 
temperature (10°C) from various tem- 
peratures from 300°C to 600°C in 
magnetic field of 0.5 Oe. are shown in 


Fig. 4. As is shown in this figure, the . 
Fig. 3. Thermo-remanent magnetization, 


direction of T.R.M. thus produced are J n-05 of therock sample. 


opposite to that of applied magnetic 
field. 
In the ordinary rock samples or ordinary ferromagnetic minerals, T.R.M.. is 
always “ positive,” ie. the direction of T.R.M. is parallel to that of applied magnetic 
She field. Such the ordinary case will ‘be 
named hereafter ‘“‘ normal thermo-rema- 
nent magnetism” or shortly ‘‘ Normal 
T.R.M.” In the case of Eig. 3, however, 
T.R.M. has a negative value, i.e. the 
direction of T.R.M. is reversed. So 
that such the case of T.R.M. will be 


Fig. 4. Thermo-remanent magnetization, Jpn-o5 called “reverse thermo-remanent 
of the separated ferromagnetic minerals, im ysis 9103 stint aatig: 


be 
* 
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magnetism ”’ or briefly ‘‘ Reverse T.R.M.”’ 


The result shown in Fig. 3 will give us a distinct and direct evidence for the 
existence of the phenomenon of self-reversal of thermo-remanent magnetism of rocks. 

On the other hand, T.R.M. of ensemble of ferromagnetic minerals separated 
from the rock is given in Fig. 4. 

In the case of Fig. 4, T.R.M.s produced by cooling from a temperature above 
300°C in magnitic field of 0.5. Oe. are always Reverse T.R.M. 

Since the content of the ferromagneic minerals in their mother rocks is less. 
than 2% in volume and those minerals grains are distributed fairly uniformly in the 
other non-ferromagnetic minerals, and further since the phenomenon of Reverse T.R.M. 
appears both in the rock samples themselves and ensemble of ferromagnetic grains 
Separated from their mother rocks, we can conclude that the Reverse T.R.M. is 
principally a magnetic character of each grain of ferromagnetic minerals. In other 
words, the self-reversal of T.R.M. in the present case ought to be due to some inner 
magnetic constitution of the ferromagnetic mineral. 

For the purpose of examining the above conclusion more directly, T.R.M. of 
special test samples, in which the contents of the concerning ferromagnetic minerals 
were diluted by adding a fair amount of magnetia were measured. The results showed 
that Reverse T.R.M.s take place in these samples regardless of the amount of content 
of the ferromagnetic minerals, so far as the temperature, from which the samples are 
cooled in a magnetic field, exceeds 300°C. 

Then, magnetic behaviours of ensembles of the ferromagnetic grains will be 
dealt with especially in detail in the following paragraphs for the purpose of examin- 
ing the characteristics of the particular phenomenon of Reverse T.R.M. 

4. Characteristics of Reverse Thermo-Remanent Magnetism of the Ferromagnetic 
Minerals. 

The relation of the partial thermo-remanent magnetism (P.T.R.M.) of the ferro- 
magnetic minerals to the temperature range for partial magnetic field cooling is 
illustrated in Fig. 5. 

The definition of partial thermo- 


remanent magnetism (6) is the follow- 
ing: During the process of cooling of 
sample from a temperature higher than 
the hightest Curie point, magnetic field 
His applied on the sample only during 
a certain temperature range from Ty+41 
to T, (where T,1>T,), the external 
magnetic field being kept zero in the 
other temperature ranges. 

The sample after thus cooled 


Fig. 5. P.T.R.M. Jin, H=2.0 (Trtt— Tn= 
‘ 50°C), of the separated ferro- 
down to the room temperature main- magnetic minerals. 


_. taines a remanent magnetism, which is ane : ; rane a0 
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called the partial thermo-remanent magnetism, and denoted by J7",, 7. In the experi- 
ment given in Fig. 5. the temverature range 4T=T,,,:—T,, was kept 59 degrees and 
the intensity of external magnetic field applied during the temperature range was 
always 2.0 Oe. 

In the ordinary rock samles or ordinary ferromagnetic minerals, J7”,, , is al- 
ways “positive.” Such the positive partial thermo-remanent magnetization will be 
called “ normal partial thermo-remanent magnetism,” or shortly normal P.T.R.M. 

In the case shown in Fig. 5, J7%,,, » is positive when T,,4:<250°C and when T,, 
>450°C, i.e. it is N.P.T.R.M. but it is distinctly negative during the temperature range 
from: T,,>250°C to T,,41<450°C. The latter case where the direction of P.T.R.M. is 
opposite to that of H will be named hereafter ‘‘ reverse partical thermo-remanent 
magnetism,’’ or shortly ‘‘ Reverse P.T.R.M.’’ 

P.T.R.M. in the magnetic field of 0.5 Oe. was also measured, the results being 
shown in Fig. 6: The general tendency of the relation’ between P.T.R.M. and the 
temperature range of magnetic field cooling in the present case is similar to that in 
the case shown in Fig. 5. 

On the other hand, it was proved experimentally that the above-mentioned relation 
of Reverse P.T.R.M. to the température range depends on the intensity of applied 
magnetic field. P.T.R.M. produced by cooling from 350°C to 250°C in various magnetic 
fields is shown in Fig. 7, where it is noticed that P.T.R.M. is reverse provided that 
the magmetic field H is less than 7.0 Oe., while it becomes normal if H exceeds 7.0 Oe. 
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5. Tempervture Coefficient of R.T.R.M. around the Ordinary Atmospheric Temperature. 
As is discussed in the later paragraph, the temperature cofficient of R.T.R.M. 


thus produced will be a significant element which is related to the constitution of this 


peculiar remanent magnetization. 


The change of the intensity of T.R.M. with temperature in a non-magnetic space 


was observed by means of a ballistic method apparatus within the temperature range 


from —70°C to 100°C. 


jax 


a. Change of J3s0,7-0,5 (7) around the 
room temperature. 
b. Change of Jgso,zr=0,5(T) around the 
room temperature. 


The results of measurements are illust- 
rated in Fig. 8, where curve (a) shows 
J20.n-02 Of the fresh ferromagnetic mineral, 
while curve (b) Ji) 7-05 Of the minearals 
which were once heated up to 800°C, the 
origina! mineral constitution being some- 
what altered. The thermal changes shown 
here are reversible at least during the 
range of temperature from —70°C to 80°C, 
and the temperature coefficient is a posi- 


tive constant in the said temperature 


renge. An interesting fact given by curve (b) is that the thermo-remanent magnetism 


which is apparently normal at the room temperature becomes reverse below —15°C. 


Generally speaking, P.T.R.M. is a function of temperature JT. Then, we can 


conclude in the above cases that 


0 
aPrJA(T)>0 


regardless of the sign of Ape 7 at the 
room temperature, so far as concern- 
ing the reversible thermal-change of 
Jz q» where the temperature change 
T.~Ty, includes the range for reverse 
P.T.R.M. shown in Fig. 5. i 

On the the other hand, Fig 9 
shows the reversible thermal changes 
of J) 7-05 and Jeo, n-0.5, Where it must 
be noticed that these two caSes were 
especially examined as the representa- 
tives of the parts of N.P.T.R.M. below 
250°C. and above 450°C respectively 
in Fig. 5 and that the magnetizations 


> x 


As is shown in the figure, the 
temperature coefficients of those 


oe © 5O> 


_magnetizations are a little. nagative, . 


— 


(1) 


Jx10° 


Fig. 9. 
Change of Jo, #-05(7) around the room temerature. 
Change of J20)' 300, 37--05(7) around the room temper- 
ature. OTE 
Change of Ji, 7-0( 7) around the room temperature. 
Change of J220: 100 95(T) around the room tempr- 
ature. ro 
Change of J), 7~-05(T) around the room temperature. 
Change of Fred pas n=05(T) around the room temper- 
ature. 
Change of Jeo’ 330’ 200, #7-:05() around the room tem- 
perature, iba weeroaas PL eee 
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but its absolute magnitude is far small compared with that shown in Fig. 5. In-other 


expressions, 
pal. 
FT ohn (PSDP RTS 
0 ot s 
oT p", H (T)20 T i> Ds (2) 


where 7, and 7,. are the critical temperatnres defined in Fig. 5 and (1). 

Since the thermal change of magnetization mentioned above is reversible with 
respect to temperature in a non-magnetic space, the magnetization concerned ought 
to be the sum of spontaneous magnetization of domains which are fixed by the process 
of magnetic field cooling. Then, we may assume that the spontaneous magnetization 
in the case of Fig. 8 contains a particular element quite different from that in the 
case of Fig. 9. In other words, it seems that the ferromagnetic mineral is com osed | 
at least of two different kinds of ferromagnetic substance. | 

In order to examine the composition and possible relations among the elemental 
magntizations the temperature coefficient of T.R.M. of the same sample treated with 
various combinations of the three representative processes of magnetic field cooling | 
was measured, where the said three processes were the cooling in magnetic field of © 
0.5 Oe. from 690°C to 500°C, from 350°C to 250°C and from 200°C to 100°C. The q 
results of experiments are shown in Table II, where, for instance, the notation J¢v} €00200,0,5 


means T.R.M. produced by the superposition of two processes of magnetic field cooling a 
from® 00°C to 500°C and from 200°C to 100°C in 0.5 Oe. 


Table II 


| Treatment of magnetic 
field cooling 


Character 


Sign of | 
of T.R.M. + 


temp. coe 


Intensity 
at el5ee 


Mean temperature coefficient 
(—70°C~+80°C) 


emu, /e om u 
+ 3.9x10-2- = RAE [degree 


| —16.7x10-2 | +105 20-4 £™"* /acgree 


Normal 


Reverse ~ 


+ 21x10% | -—2 5107 SB ame Normal | 


= F00;2005 2-0,5 + 2.1107 — 2 x10-6 SU" /aegree 


i ee —20.7x10-2 | 410. 0x 10-4 £m faegree 
—T41x10% |. + 9.6X 10-1 2" /aesree 


—16.3X 10-2 


ih 10.5% 10-4 25 /aegsee 28 
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(c) The addition law of P.T.R.M., which has been established for the case of 
ordinary T.R.M., does not exactly hold in the present case. 
6. Process of Development of Reverse Thermo-Remanent Magnetism. 

In order to see directly the process of development of R.T.R.M., the change in 
: magnetization of the ferromag- 
netic minerals with temperature 
was measured by means of a 
ballistic method schematically 
searck coil illustrated in Fig. 10. The 
results of experiment are given 


Ss. 


_ Fig. 10. Measuring apparatus for development and , 
decay of T.R.M. in Figures 11~19 where the hollw 
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Fig. 16. Development and decay of Jz, pry cee 
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addition of the reversible magnetization /,. 


In the re-heating process without magnetic field, on the contrary, the measured 
magnetization is to be shown by 


J(I)}= _% Zo nt eT) the 9 8 
_ al T) LTS Eas 


In the series of Figures 11~19, the changes with temperature of J2%(T), J#(T), 
Ji0(T), F(T), FR(T), FET), PRT), J#&(T) and Ji (T) are illustrated successively. 

The suffix concerning H in the notation of Je w(T) will be omitted hereafter, 
since the magnitude of H is kept constant throughout all experiments dealt with in 
this paragraph, i.e. H=0.5 Oersteds. 

The magnitude of Ir i A(T) in the cooling process and that in the re-heating 
process given in these Roukes show a definite amount of discrepancy. We may assume 
at the present stage that the said discrepancy is due to the condition that the mag- 
netization in the cooling process contains a little after-effect of foregoing magnetization 
with respect to time. In other words, it seems that the magnetic domains were not 
perfectly fixed a little while after removal of external magnetic field in the cooling 
process. 

J&(T), J2(T) and J%(T) belong to Ni-process described in § 5, and they are 
normal in the direction of magnetization. However, J3x(T) and Ji(T) rise up around 
450°C both in the cooling and re-heating processes. This is a particular phenomenon, 
since » Fe (Zl) curve of the ordinary T.R.M.- increases monotonously with decrease in 
temperature in accordance with the Curie’s law. 

“0(T) in Fig. 12 shows a more complicated form, the remanent magnetization 
decreasing with decrease in temperature below 200°C. ‘This is also a particular case 
compared with the characteristics of the ordinary T.R.M. 

J2\(T) shown in Fig. 14 is a transient state between Ni-process and &-process, 
giving very little remanent magnetization. 

On the other hand, J20(T), J3(T) and J3%(T) belong to R-process mentioned in 
§ 5. The most typical example of this group will be /3(T) shown in Fig. 16, where 
the remanent magnetization decreases almost linearly with decrease in temperature 
after the removal of external magnetic field at 250°C, and crosses the zero-line at 
about 200°C, becoming the reverse T.R.M. below the said critical temperature. The 
change in J¥*(T) with temperature is almost-reversible so far as T’ does not exceed 

200°C. It must be recalled that the above-mentioned reversible and linear change with 
temperature holds in the temperature range, at least, down to —70°C. 

The general tendency of the curves of J%(T) and Jz#(7) are the same as that 

of J%(T), though the reverse T.R.M. in the former two cases are small. . 
Je(T) and J%(T) given in Figures 18 and 19 are generally of the normal form, 


-and these cases are classified into N.-process in § 5. 
4 
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It may be concluded from the above-mentioned results that the reverse T.R.M. 
is produced chiefly in the temperature range from 300°C to 250°C or in somewhat 
wider range under the effect of an external magnetic field. 

Further, it must be noticed by summarizing the results given in § 5 and § 6, 
that the positive temperature coefficient of the reverse T.R.M. is fairly large and kept 
nearly constant through a fairly wide range of temperature, since this is to be a 
distinct contrast against the negative temperature coefficient of ordinary T.R.M. which 
follows the Curie’s law. 

It will also be worthwhile to point out that the rise up of J A (T) curve given 
by M,-process in a certain temperature range can hardly harmonize with the ordinary 
conception of ferromagnetism, so far as we assume that T.R.M. is the sum of spont- 
aneous magnetizations of the domains fixed by the process of magnetic field cooling. 


7. Process of Development of Total Thermo-Remanent Magnetization in the Case where 
Reverse P.T.R.M. is enclosed. 


In the view of geophysical standpoint, one’s interest may concern the process of 
development of reversed total thermo-rement magnetization which is produced by — 
cooling from a certain high temperature throughly in a magnetic field. Hence, the — 
26 Ixlo'emuZ, a5 above-mentioned process was meas- , 
aw. ured by the apparatus described 
in the preceding paragraph. Fig. 20 q 
shows an example of measure- — 


ments, where H was kept 0.5 Oc. a 


throughout the whole process of 
cooling. In this figure, the change 

if 
in magnetization with temperate @ 
is also shown in the. case of initial E 
heating of the same “sample in e 
same magnetic fel where the 


beans 5. SS Sr 
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between the J.(T) curve and the J,(T) curve at a same temperature in Fig. 20 is shown 
by 


FAT IAT VA TEAT) (8) 


The difference J,(T)—J,(T) thus derived is also shown in Fig.20 by the borken line 
curve which becomes negative below 85°C. 

After the removal of H at room temperature (10°C), the sample has a negative 
remanent magnetization, i.e. Reverse T.R.M., as already described in ¢ 3. And the 
intensity of the Reverse T.R.M. is just coincident with the magnitude of J,(10°C)—/, 
(10°C). 

The sample thus treated was again heated in a non-magnetic space, the change 
of T.R.M. in this case being shown by full circles in Fig.20. The magnetization in 
this case is to be equal to of ja(l), provided that T.R.M. is assumed to be caused by 
fixing a certain part of magnetic domains by the magnetic field cooling. It will be 
seen in the figure that JZ, (T) thus directly measured is almost equal to J,(T)—J,(T) 
at any temperature. This fact may indicate that the relation given by (8) and con- 
sequently that given by (7) holds fairly well in the present case. 

8. Possible Mechanism of Causation of Reverse Thermo-Remanent Magnetization. 

It will be concluded from the experimental results that Reverse T.R.M. is a 
definite physical phenomenon, reproducible and universal, provided that certain physical 
conditions are fulfilled in the ferromognetic grains. There will be a few theoretical 
possibilities of self-reversal of thermo-remanent magnetization, as already been pointed 
out by L. Néel (10). Those possibilities may belong to two different categories. 

The first one is that the ferromagnetic mineral sample contains more than two 
different constituents which have an intense magnetic interaction to each other, so 
that one of them is magnetized by the demagnetizing field of another which has been 
magnetized by the external field, the remanent magnetization of the former bringing 
about Reverse T.R.M. according as the specimen is cooled. 

The other case is the ferrimagnetic substance, spontaneous magnetization of 
which belongs to N- or V-type according to Néel’s classification (13). . 

It does not seem that our experimental results obtained so far are enough 
sufficient to allow one to determine finally which one of these mechanisms should be 
the real. The following discussions are to treat this problem briefly on the basis of 
the experimental results, where are used the notations that . 

T,....domain fixing temperature, 
H.....intensity of external magnetic field, 


J, ....irreversible or fixed magnetization, 
J, ....reversible magnetization, AG 
o ....intensity of spontaneous magnetization, 


n....number of fixed magnetic domains per unit volume, i 
. v ....volume of a magnetic domain. 
The capital letters such as A, B, and C represent the ferromagnetic constituents 


ee 
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involved. 
One may notice from Figs. 11, 12, 13, 18 and 19, that the magnetizations produced 


in N; and N, processes are not playing any essential role in the formation of Reverse 
T.R.M., although the abrupt increase of magnetization at about 450°C in Figs. 12 and 
20 will be an interesting fact. Therefore the domain fixing temperatures of the 
constituents essential to Reverse T.R.M. must be in the comparatively narrow range, 
i.e. the range where the FR-process takes place. 
1) Reverse T.R.M. as a result of magnetic interaction between different constituents. 
The first type of the theoretical possibilities is based upon following assumptions. 
The sample contains two different ferromagnetics A and B having different domain 
fixing temperatures 7,4 and 7,” and the magnetic domains of A are fixed at 7,“ in 
the direction of H during the process of cooling in H, resulting in an irreversible 


magnetization as follows. 


SJ AMTI=n 0, OT), “where © TXT,A. | (9) 
Then the effective field H.yy on the domains of B at T=7T,? becomes 
Heyy = H—-{D-JAT,?)+a-JA(T.7)+8-J7(T.7)} (10) 


where D, a, and # are some constants representing demagnetizing factors. If TAT) 
increases as much as to make His in (19) negative, the total irreversible magnetization 
at T<T,” becomes = 

FAT)H=T AT) AT AD) 2104-04-61) p- 0g (T) (11) 
Hence, if the condition ; 


Nava — 9,(T) 
AyD p sae be 


; 


holds, the resultant magnetization can be negative. Moreover, since I(T) in (11) EG q 


theoretically or experimentally, positive at hg the condition to make it negative at 
room As becomes 


rae . bint oT) a TY ic bosehe A od Te ake ae 
' OAT 8) < Ox( Ly") * ferns iret Bae 3) : 


This Perditon seems . plausible from our panne sale that ae temperature 
coefficient of J,(T’) produced in R-process is large. Thus, if o,(T) increase 
gee IAT) at the room. temaerstvinnnts sade TM 06 eee 


= 
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where d and s denote the difference in demagnetizing factors of sphere and the ellip- 
soid in the direction of long axis and the susceptibility as to the irreversible magneti- 
zation at T,. In (14), the first and second terms correspond to the irreversible 


magnetization of B and A, and condition for appearance of Reverse T.R.M. is, 


Trt eid y<U, T....room temperature. (15) 


Néel derived a more explicit formula from (14) by using a simpler model and shown 
that the occurence of Reverse T.R.M. is not impossible provided that proper values 
are assumed for the quantities involved in the equation. 

The above mentioned Néel’s model seems to be consistent in itself and actually 
probable. This model, however, does not seem to be applied for the explanation of 
our experimental results: if the above theory is accepted, P.T.R.M., de zi fn siwhere 
T1<T,?<T:<T,4, must be in the normal direction and have a large temperature 
coefficient which is minus in sign, because in this model the observed large” positive 
coefficient of Reverse T.R.M. is to be caused by the steep development of o;(T) with 
decrease in temperature. But no such the result has ever been obtained. Another 
example unfavorable to this machanism will be that the re-heating curves of irreversible 
magnetization produced by R-process hardly exceed the zero line, for this tendency 
makes it quite dubious that there were such an intense normal remanent. magnetization 
at 7,” as to magnetize B in the reverse direction by its demagnetizing field in the 
cooling process. 

In the present paper, we shall tentatively attack another possibility. The above 
facts and another that Reverse T.R.M. is fairly stable, changing reversibly with 
temperature below 7:, may suggest that the whole process of fixing of domains in 
reverse direction takes place during the cooling from T, and T;. If so, the domains 
in question must be orientated into the reverse direction while the extertal magnetic 
field magnetizes the total magnetization, including the reversible part of it, in the 
normal direction. Examing /(T)-7T curves in Figs. 16 and 17, for example, we can 
observe the fact that the reversible part of magnetization is markedly large. This 
means that the reversible part should be the major factor in producing the demagnet- 
izing field. 

Now we shall assume a somewhat modified model, where the sample is composed 
also of constituents A and B. Constituent A is the main constituent of this sample 
and is supposed to be the same with the normal ferromagnetic minerals. o(7) in 
Fig. 2 would be the spontaneous magnetization of A. The magnetization of A is 
composed of two parts as is in the other ordinary specimens, i.e. the reversible part 7 
a, the irreversible part C. The domain fixing temperature of constituent B is to be 
in the R-temperature range and o,7(T) incerases so steeply as to cause the large 
temperature coefficient of Tet n(T) observed. B is assumed to be ina peculiarly close 
magnetic coupling with 7. The content of B is scanty compared with that of A. 
Therefore J,” is easily saturated while J, and J;° increase with the increase of H. In 


the remanent magnetization by stronger field, Jif begins to play the major role, 
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According to the model described above, let us investigate Tos 7 produced in 
R-process, ie. T:<T,?7<T:. Magnetization is expressed in general as, 
KT )=JA(T)+I (1). (16) 
When H is weak, the equation (16) becomes, 
ID) =F ADT) tT) =2(L):o(T)0-H+Ji(T), where T,?<T<T; (17) 
KT)=JAT) AI AT) tI AT) 
=n,(T)-6,(T)+0,-H+22(Ty*)-09(T) 0x A— Der TAT» ASAT) 


where Lynd <a; (18) 
KT) =J.2AT) 4+ F(T) = 2 F(Ty?)-0(T)-08{H— Dor SAT s®)} + IAT) 
where T<T\. (19) 


In these equations 2 denotes the number of magnetic domains, in a unit volume, 
susceptible to a unit magnetic field and D;, represents the demagnetizing factor of 
Jeu a affecting on B. If one puts J,(7)=0 in the above equations, the expressions for 
the PeResttie curves are to be obtained. Hence, as long as J,°(7) is small, the tend- 
ency that re-heating curves do not exceed the zero line can be explained, and the 
necessary condition for Reverse T.R.M. is > 

H— Dar FAT y") <0. (20) 

In this model, the essential problem is how such an intense magnetic interaction — 
can be possible. The macroscopic model as Néel’s is not applicable in our case, because — 
‘the concentration of ferromagnetic grains does not affect the magnitude of Reverse 4 
T.R.M. Such a close magnetic interaction should be, in the authors’ opinion, of the | 


domain scale. 
When dH is strong, the magnitude of I(T) increases and the reheating curves — 
begin to exceed the zero line and at last J,;(T) at room temperature becomes normal. 
2) A model containing N- or V-type ferrimagnetic substance. 
pees has been mentioned in the beginning of this paragraph, another almost 
equally probable mechanism may be ‘establis! Bae 
if the existence of N- or V-type ferrimag 
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negative below 9,. This model, however, is obliged to have the same difficulty with 
the model containing two fixing constituents in explaining the observed distinct tendency 
that the re-heating curves in Figs. 16 and 17 hardly exceed the zero line. 

9. Conclusion and Further Problems. 

It seems that the existence of self-reversal of thermo-remanent magnetism in 
nature was proved by the present study, and that Reverse T.R.M. can be attributed 
to the magnetic characteristics of ferromagnetic (or ferrimagnetic) mineral grains. 

Physical mechanism of causation of this particular phenomenon, however, has 
not yet been fully clarified, though several facts directly connected to the mechanism 
were obtained in the present experiments. The problem of Reverse T.R.M. will be 
significant from two different standpoints of view, one being in connection with geo- 
physical problem, especially with palaeomagnetism, while the other is with physics of 
magnetism, especially with characteristics of ferrimagnetic substances. 

From the standpoint of geophysics, it will be desirable that a large number of 
igneous rocks which have reverse remanent magnetization zm situ are examined in the 
same way as here. Studies of some of other Japanese samples are now under going. 
But it will be further desirable that the rook samples collected from various localities 
over the world are examined by a systematic method. 

From the standpoint of physical study of “ferrite” substances, the mechanism 
of development of Reverse T.R.M. will be especially interesting, since the phenomenon 
looks like closely relating to Néel’s ferrimagnetism. The best way will be to separate 
out the component of the mineral grains which are directly responsible for Reverse 
T.R.M., and to study various physical characteristics of that compdnent. Another 
problem will be the behavior of Reverse T.R.M. at lower temperatures. The study of 
the latter problem is now under working. 

In concluding, the writers wish to express their sincere thanks to Dr. H. Kuno 
and Dr. I. Iwasaki for their kind cooperation in petrological and chemical works res- 
pectively in the present study. Their hearty thanks are also due to Prof. L. Néel for 
his valuable dicussions in the course of this study, and to Prof. S. Nakamura and Prof. 


M. Hasegawa for their interest and encouragement. 
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LETTER TO THE EDITORS 


Semi-Diurnal Lunar Variations in the Sporadic-E 


The critical penetrating frequency of the sporadic-F, fEs, observed at Kokubunji 
(139°29.3’E, 35°42.4’N) was analyzed for 12 months of June solstices in 1949-51 to obtain 
the semi-diurnal lunar variations in the sporadic-E. For the purpose of eliminating 
the’solar variations, difference (4fEs) between hourly value of fEs and corresponding 
monthly median value was calculated. And, they were arranged in lunar time, reckoned 
from upper transit, on one sheet during 29 days [1] after the epoch of the full and 
the new moon, and their mean values were obtained. Same treatments were also 
conducted for 2004fEs/medianfEs [2]. 

Fig. 1 is the mean 
diurnal curves of 9 mean 
lunar variations for 4fEs 
and 200 4fEs/median fEs 
obtained from these treat- 
ments. The semi-diurnal 
variations will be obvious - 
in the figure. 


A more accurate 
study is provided’ by 


harmonic analysis of 
these variations. Reliable 


- 1p ; : = ‘aa nine results, by the analy- 


Fig. 1 Average lunar variation curves of AfEs (thick full- sis of the statistical signi- 
line) and 200 AfEs/median /Es (thick dotted-line) at ficance [3], are shown by 
Kokubunji during summer months for the year 1949- 
51. Almost all errors are contained in the band 
between the two thin lines. (confidence limit,95.5%) (Both G- and E-symbols 


in fEs data are assumed 

to be 0 and 0.9 MC/S, and mean values of the results for these two cases are plotted 

in the figure.) On the average, the lunar semi-diurnal variation attains its maximum 

amplitude of about 0.35 MC/S and of about 13% for 2004fEs/medianfEs, approxi- 
mately 8 hours after the lunar culminations. z. 

The amplitude of the diurnal variation is generally smaller than the semi-diurnal, 

and-the phase distributes in all quadrants of the harmonic dial. The lunar variation: 


harmonic dials in Fig. 2. 


in December solstice and equinoxes is difficult to obtain, because of few occurrence 
of Es. 
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Fig. 2 Harmonic dials of semi-diurnal lunar varia- 


tion in the Es for AfEs (left) and 200AfEs/ 
medianfEs (right) at Kokubunji during 
summer months for the year 1949-51. 
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Though the individual values of 
fEs data may have not enough 
accuracy to enable one to draw any 
tidal oscillation, the results of the 
above statistical study might be 
considered as giving some informa- 
tion on the semi-diurnal lunar tide 
in the Es region. 

It will be attractive that the 


amplitude of the lunar 12 hours. 


variation in the Es is much larger 


than that in the & and even in the” 
F2 [4], and that the phase is resem- | 


ble to the F2, though it is different 
from the phase of the lunar 
variation of the absorption in the 
D region [5]. 

If these results might be consi: 


dered in the light of Martyn’s 


eds 4 


(i: 1949, 4: 1950, @:1951, x : mean) theory [6] of drift movements with 


horizontal tidal winds, a suggestion 
for the formation of the Es might be found. In conclusion, the writer wishes to express y 

his sincere thanks to Prof. M. Hasegawa for helpful suggestions, and also to Miss 
M. Kato and Miss T. Nakamura for their assistances. (Read: April 28, 1952) 
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